A calf with hind limb paralysis and dysstasia and a genome sequence analysis of an isolated Clostridium perfringens toxinotype E strain.

Clostridium perfringens toxinotype E infections are rare in calves, and the development of intestinal lesions were commonly observed. In 2012, a 6-day-old calf in Japan exhibited swelling with emphysema on the right gluteal region, sudden paralysis of the hind limb and dysstasia. A pathological examination revealed myositis of the gluteal muscle and neuritis of the ischiatic nerve. C. perfringens type E strain CP118 was isolated from the affected muscle. However, the intestinal symptoms and lesions that commonly develop in type E infections in calves were not detected in the present case. Genome analyses revealed that CP118 possessed 16 virulence-related genes, including enterotoxin, and was closely related to other type E and F strains. Particularly, CP118 was more closely related to type E strains from humans, including a food poisoning case, than calf isolates, suggesting its potential to cause food poisoning in humans and, thus, its importance as a potential risk to public health. Since CP118 did not possess the reported toxin genes associated with neuropathy, pyogenic inflammation caused by CP118 and/or other bacteria may have damaged the ischiatic nerve, resulting in neuropathy. Alternatively, unidentified CP118 toxins may have caused the neuropathy. This is the first study to report C. perfringens type E infection with peripheral neuropathy. The distribution of all the reported virulence-related genes in the C. perfringens population as well as the details of this rare case will provide further insights into C. perfringens type E infections.

Fine mapping of Rf5 region for a sorghum fertility restorer gene and microsynteny analysis across grass species.

Cytoplasmic male sterility (CMS) is widely used to control pollination in the production of commercial F1 hybrid seed in sorghum. So far, 6 major fertility restorer genes, Rf1 to Rf6, have been reported in sorghum. Here, we fine-mapped the Rf5 locus on sorghum chromosome 5 using descendant populations of a 'Nakei MS-3A' × 'JN43' cross. The Rf5 locus was narrowed to a 140-kb region in BTx623 genome (161-kb in JN43) with 16 predicted genes, including 6 homologous to the rice fertility restorer Rf1 (PPR.1 to PPR.6). These 6 homologs have tandem pentatricopeptide repeat (PPR) motifs. Many Rf genes encode PPR proteins, which bind RNA transcripts and modulate gene expression at the RNA level. No PPR genes were detected at the Rf5 locus on the corresponding homologous chromosome of rice, foxtail millet, or maize, so this gene cluster may have originated by chromosome translocation and duplication after the divergence of sorghum from these species. Comparison of the sequences of these genes between fertile and CMS lines identified PPR.4 as the most plausible candidate gene for Rf5.

A feasibility trial of genomics-based diagnosis detecting insecticide resistance of the diamondback moth.

BACKGROUND: Insecticide resistance management has been key for crop protection for over 70 years and is increasingly important because the development of new active ingredients has decreased in recent years. By monitoring the development of resistance in a timely manner, we can effectively prolong insecticide efficacy. Genomic-based diagnosis can reliably predict resistance development if information on resistant mutations against major pesticides is available. Here, we developed a feasibility trial of genomics-based diagnosis of insecticide resistance in diamondback moth (Plutella xylostella) populations in Nagano Prefecture, Japan. Amplicon sequencing analyses using a next-generation sequencer (Illumina MiSeq) for major insecticides, including diamides, pyrethroids, Bacillus thuringiensis (Bt) toxin (Cry1Ac), organophosphates, and spinosyns, were conducted. RESULTS: Mutations related to the resistance of pyrethroids, organophosphates, and diamides (flubendiamide and chlorantraniliprole) prevailed, while those of a diamide (cyantraniliprole), Bt (Cry1Ac), and spinosyns were scanty, suggesting that they are still effective. The results of the genomics-based diagnosis were generally concordant with the results of bioassays. Resistance development tendencies were generally uniform across Nagano. CONCLUSION: An insecticide-resistance management campaign can be conducted in Nagano Prefecture with a quick genomic-based diagnosis in early spring while bioassay is the only option for monitoring resistances whose mutations are unavailable. Our study is the first step in the future management of insecticide resistance in all significant pests. © 2022 Society of Chemical Industry.

Development of an SSR marker-based genetic linkage map and identification of a QTL associated with flowering time in Eustoma.

Lisianthus (Eustoma grandiflorum) is an important floricultural crop cultivated worldwide. Despite its commercial importance, few DNA markers are available for molecular genetic research. In this study, we constructed a genetic linkage map and to detect quantitative trait loci (QTLs) for important agronomic traits of lisianthus. To develop simple sequence repeat (SSR) markers, we used 454-pyrosequencing technology to obtain genomic shotgun sequences and subsequently identified 8263 putative SSRs. A total of 3990 primer pairs were designed in silico and 1189 unique primer pairs were extracted through a BLAST search. Amplification was successful for more than 1000 primer pairs, and ultimately 278 SSR markers exhibited polymorphism between the two lisianthus accessions evaluated. Based on these markers, a genetic linkage map was constructed using a breeding population derived from crosses between the two accessions, for which flowering time differed (>140 days when grown under 20°C). We detected one QTL associated with flowering time (phenotypic variance, 27%; LOD value, 3.7). The SSR marker located at this QTL may account for variation in flowering time among accessions (i.e., three accessions whose nodes of the first flower were over 30 had late-flowering alleles of this QTL).

Multiple Wheat Genomes Reveal Novel Gli-2 Sublocus Location and Variation of Celiac Disease Epitopes in Duplicated α-Gliadin Genes.

The seed protein α-gliadin is a major component of wheat flour and causes gluten-related diseases. However, due to the complexity of this multigene family with a genome structure composed of dozens of copies derived from tandem and genome duplications, little was known about the variation between accessions, and thus little effort has been made to explicitly target α-gliadin for bread wheat breeding. Here, we analyzed genomic variation in α-gliadins across 11 recently published chromosome-scale assemblies of hexaploid wheat, with validation using long-read data. We unexpectedly found that the Gli-B2 locus is not a single contiguous locus but is composed of two subloci, suggesting the possibility of recombination between the two during breeding. We confirmed that the number of immunogenic epitopes among 11 accessions varied. The D subgenome of a European spelt line also contained epitopes, in agreement with its hybridization history. Evolutionary analysis identified amino acid sites under diversifying selection, suggesting their functional importance. The analysis opens the way for improved grain quality and safety through wheat breeding.

Antimicrobial Resistance Genes in Bacteria Isolated From Japanese Honey, and Their Potential for Conferring Macrolide and Lincosamide Resistance in the American Foulbrood Pathogen Paenibacillus larvae.

American foulbrood (AFB) is the most serious bacterial disease of honey bee brood. Spores of the causative agent Paenibacillus larvae are ingested by bee larvae via brood foods and germinated cells proliferate in the larval midgut. In Japan, a macrolide antibiotic, tylosin, is used as the approved prophylactic for AFB. Although tylosin-resistant P. larvae has yet to be found in Japan, it may emerge in the future through the acquisition of macrolide resistance genes from other bacteria, and bacteria latent in brood foods, such as honey, may serve as a source of resistance genes. In this study, to investigate macrolide resistance genes in honey, we attempted to isolate tylosin-resistant bacteria from 53 Japanese honey samples and obtained 209 isolates from 48 samples in the presence of 1 µg/ml of tylosin. All isolates were Gram-positive spore-forming bacteria mainly belonging to genera Bacillus and Paenibacillus, and 94.3% exhibited lower susceptibility to tylosin than Japanese P. larvae isolates. Genome analysis of 50 representative isolates revealed the presence of putative macrolide resistance genes in the isolates, and some of them were located on mobile genetic elements (MGEs). Among the genes on MGEs, ermC on the putative mobilizable plasmid pJ18TS1mac of Oceanobacillus strain J18TS1 conferred tylosin and lincomycin resistance to P. larvae after introducing the cloned gene using the expression vector. Moreover, pJ18TS1mac was retained in the P. larvae population for a long period even under non-selective conditions. This suggests that bacteria in honey is a source of genes for conferring tylosin resistance to P. larvae; therefore, monitoring of bacteria in honey may be helpful to predict the emergence of tylosin-resistant P. larvae and prevent the selection of resistant strains.

Investigation of the Genetic Diversity of a Rice Core Collection of Japanese Landraces using Whole-Genome Sequencing.

The Rice Core Collection of Japanese Landraces (JRC) consisting of 50 accessions was developed by the genebank at the National Agriculture and Food Research Organization (NARO) in 2008. As a Japanese landrace core collection, the JRC has been used for many research projects, including screening for different phenotypes and allele mining for target genes. To understand the genetic diversity of Japanese Landraces, we performed whole-genome resequencing of these 50 accessions and obtained a total of 2,145,095 single nucleotide polymorphism (SNPs) and 317,832 insertion-deletions (indels) by mapping against the Oryza sativa ssp. japonica Nipponbare genome. A JRC phylogenetic tree based on 1,394 representative SNPs showed that JRC accessions were divided into two major groups and one small group. We used the multiple genome browser, TASUKE+, to examine the haplotypes of flowering genes and detected new mutations in these genes. Finally, we performed genome-wide association studies (GWAS) for agronomical traits using the JRC and another core collection, the World Rice Core Collection (WRC), comprising 69 accessions also provided by the NARO genebank. In leaf blade width, a strong peak close to NAL1, a key gene for the regulation of leaf width, and, in heading date, a peak near HESO1 involved in flowering regulation were observed in GWAS using the JRC. They were also detected in GWAS using the combined JRC + WRC. Thus, JRC and JRC + WRC are suitable populations for GWAS of particular traits.

Origin, selection, and spread of diamide insecticide resistance allele in field populations of diamondback moth in east and southeast Asia.

BACKGROUND: The investigation of molecular mechanisms and evolution of resistance to insecticides is an ongoing challenge, as researchers must provide guidance to manage the resistance to achieve sustainable production in agriculture. Predicting, monitoring, and managing insecticide resistance requires information on the origins, selection, and spread of resistance genes. The resistance of Plutella xylostella (L.) against diamide insecticides is becoming an increasingly severe problem in east and southeast Asia. In this study, the evolution of resistance was investigated using a resistance allele [ryanodine receptor (RyR); G4946E mutation] and its flanking regions, as well as mitochondrial cytochrome c oxidase subunit I (mtCOI). RESULTS: The sequences of the flanking region of the G4946E and mtCOI suggested that the G4946E mutation has a key role in resistance. Furthermore, the G4946E mutation has multiple origins, and congenic resistant mutations have spread across east and southeast Asia, despite substantial geographical barriers. In addition, the susceptibility of field populations partially recovered during winter, based on the observed decrease in the G4946E (resistant allele) frequency. Finally, the resistance level indexed by the frequency of the E4946 allele was significantly lower in non-overwintering regions than in overwintering regions. CONCLUSION: The information of the present study is useful to monitor resistance using molecular markers and to develop strategies to delay the evolution of diamide resistance.

Multiple wheat genomes reveal global variation in modern breeding.

Advances in genomics have expedited the improvement of several agriculturally important crops but similar efforts in wheat (Triticum spp.) have been more challenging. This is largely owing to the size and complexity of the wheat genome1, and the lack of genome-assembly data for multiple wheat lines2,3. Here we generated ten chromosome pseudomolecule and five scaffold assemblies of hexaploid wheat to explore the genomic diversity among wheat lines from global breeding programs. Comparative analysis revealed extensive structural rearrangements, introgressions from wild relatives and differences in gene content resulting from complex breeding histories aimed at improving adaptation to diverse environments, grain yield and quality, and resistance to stresses4,5. We provide examples outlining the utility of these genomes, including a detailed multi-genome-derived nucleotide-binding leucine-rich repeat protein repertoire involved in disease resistance and the characterization of Sm16, a gene associated with insect resistance. These genome assemblies will provide a basis for functional gene discovery and breeding to deliver the next generation of modern wheat cultivars.

Draft Genome Sequence of Brucella ceti Isolated in the Western Pacific Ocean.

In 2018, Brucella ceti was isolated from a bottlenose dolphin from the western Pacific Ocean. Here, we report a draft genome sequence of the isolate BD1442 of sequence type 27, which is the only sequence type known to have been isolated from human clinical cases.

Evolutionary dynamics and impacts of chromosome regions carrying R-gene clusters in rice.

To elucidate R-gene evolution, we compared the genomic compositions and structures of chromosome regions carrying R-gene clusters among cultivated and wild rice species. Map-based sequencing and gene annotation of orthologous genomic regions (1.2 to 1.9 Mb) close to the terminal end of the long arm of rice chromosome 11 revealed R-gene clusters within six cultivated and ancestral wild rice accessions. NBS-LRR R-genes were much more abundant in Asian cultivated rice (O. sativa L.) than in its ancestors, indicating that homologs of functional genes involved in the same pathway likely increase in number because of tandem duplication of chromosomal segments and were selected during cultivation. Phylogenetic analysis using amino acid sequences indicated that homologs of paired Pikm1-Pikm2 (NBS-LRR) genes conferring rice-blast resistance were likely conserved among all cultivated and wild rice species we examined, and the homolog of Xa3/Xa26 (LRR-RLK) conferring bacterial blight resistance was lacking only in Kasalath.

Domain Unknown Function DUF1668-Containing Genes in Multiple Lineages Are Responsible for F1 Pollen Sterility in Rice.

Postzygotic reproductive isolation maintains species integrity and uniformity and contributes to speciation by restricting the free gene flow between divergent species. In this study we identify causal genes of two Mendelian factors S22A and S22B on rice chromosome 2 inducing F1 pollen sterility in hybrids between Oryza sativa japonica-type cultivar Taichung 65 (T65) and a wild relative of rice species Oryza glumaepatula. The causal gene of S22B in T65 encodes a protein containing DUF1668 and gametophytically expressed in the anthers, designated S22B_j. The O. glumaepatula allele S22B-g, allelic to S22B_j, possesses three non-synonymous substitutions and a 2-bp deletion, leading to a frameshifted translation at the S22B C-terminal region. Transcription level of S22B-j and/or S22B_g did not solely determine the fertility of pollen grains by genotypes at S22B. Western blotting of S22B found that one major band with approximately 46 kDa appeared only at the mature stage and was reduced on semi-sterile heterozygotes at S22B, implying that the 46 kDa band may associated in hybrid sterility. In addition, causal genes of S22A in T65 were found to be S22A_j1 and S22A_j3 encoding DUF1668-containing protein. The allele of a wild rice species Oryza meridionalis Ng at S22B, designated S22B_m, is a loss-of-function allele probably due to large deletion of the gene lacking DUF1668 domain and evolved from the different lineage of O. glumaepatula. Phylogenetic analysis of DUF1668 suggested that many gene duplications occurred before the divergence of current crops in Poaceae, and loss-of-function mutations of DUF1668-containing genes represent the candidate causal genetic events contributing to hybrid incompatibilities. The duplicated DUF1668-domain gene may provide genetic potential to induce hybrid incompatibility by consequent mutations after divergence.

Metagenomic profiles of the rumen microbiota during the transition period in low-yield and high-yield dairy cows.

We investigated potential relationships between rumen microbiota and milk production in dairy cows during the transition period. Twelve dairy cows were divided into a low-yield (LY) or high-yield (HY) group based on their milk yield. Rumen samples were taken from dairy cows at 3 weeks before parturition, and at 4, 8, and 12 weeks after parturition. 16S rDNA-based metagenomic analysis showed that diversities of rumen microbiota in both groups were similar and the number of operational taxonomic units (OTUs) was lower in the postpartum than prepartum period in both groups. The abundance of Bacteroidetes and ratio of Bacteroidetes:Firmicutes was higher in the HY than the LY group. OTUs assigned to Prevotella bryantii, Fibrobacter succinogenes, Ruminococcus albus, Butyrivibrio fibrisolvens, and Succinivibrio sp. were abundant in the HY group. These OTUs were significantly related to the propionate molar proportion of rumen fluids in the HY group. OTUs assigned to Lachnospiraceae, Bifidobacterium sp. and Saccharofermentans were dominant in the LY group. Predictive functional profiling revealed that abundance of gene families involved in amino acid and vitamin metabolism was higher in the HY than the LY group. These results suggest that the community structure and fermentation products of rumen microbiota could be associated with milk production of dairy cows.

Scanning RNA-Seq and RAD-Seq approach to develop SNP markers closely linked to MALE STERILITY 1 (MS1) in Cryptomeria japonica D. Don.

Cryptomeria japonica is a major forestry tree species in Japan. Male sterility of the species is caused by a recessive gene, which shows dysfunction of pollen development and results in no dispersed pollen. Because the pollen of C. japonica induces pollinosis, breeding of pollen-free C. japonica is desired. In this study, single nucleotide polymorphism (SNP) markers located at 1.78 and 0.58 cM to a male sterility locus (MS1) were identified from an analysis of RNA-Seq and RAD-Seq, respectively. SNPs closely linked to MS1 were first scanned by a method similar to MutMap, where a type of index was calculated to measure the strength of the linkage between a marker sequence and MS1. Linkage analysis of selected SNP markers confirmed a higher efficiency of the current method to construct a partial map around MS1. Allele-specific PCR primer pair for the most closely linked SNP with MS1 was developed as a codominant marker, and visualization of the PCR products on an agarose gel enabled rapid screening of male sterile C. japonica. The allele-specific primers developed in this study would be useful for establishing the selection of male sterile C. japonica.

Identification of candidate genes responsible for the susceptibility of apple (Malus × domestica Borkh.) to Alternaria blotch.

BACKGROUND: The mechanism underlying the interaction between host plant and host-selective toxin (HST)-producing Alternaria alternata during infection is of particular interest for sustainable crop production. Alternaria blotch of apple (Malus × domestica Borkh.) caused by A. alternata apple pathotype is a major disease particularly in East Asia, which is the largest producer of apples globally. A single dominant gene, Alt, controls the susceptibility of the apple cultivar 'Delicious' to Alternaria blotch. In this study, we fine mapped the Alt locus and characterized three potential candidate genes. RESULTS: We used 797 F1 individuals derived from 15 crosses between apple accessions susceptible (Alt/alt) and resistant (alt/alt) to Alternaria blotch to construct physical and genetic maps of the Alt locus located on the top of chromosome 11. Susceptible accessions were derived from 'Delicious.' To fine map the Alt locus, we constructed a BAC library of 'Starking Delicious,' a sport of 'Delicious,' and used graphical genotyping to delimit the Alt locus to a region of 43 kb. Three genes predicted within the candidate Alt region were potentially involved in plant defense response, among which the gene encoding a coiled coil-nucleotide binding-leucine rich repeat (CC-NB-LRR) type disease resistance protein was the most promising. Moreover, a 12-bp insertion was uniquely identified in the 5' untranslated region of the Alt-associated allele of this gene, the presence or absence of which co-segregated with the susceptibility or resistance to A. alternata apple pathotype, respectively, among 43 tested cultivars including old ones and founders of modern apple breeding. CONCLUSION: A disease resistance protein has been suggested as a determinant of susceptibility/resistance to HST-producing A. alternata for the first time. Our finding provides new insight into the mechanism of HST-mediated disease control used by A. alternata against host plants.

Structural features of two major nucleolar organizer regions (NORs), Nor-B1 and Nor-B2, and chromosome-specific rRNA gene expression in wheat.

The reference genome sequence of wheat 'Chinese Spring' (CS) is now available (IWGSC RefSeq v1.0), but the core sequences defining the nucleolar organizer regions (NORs) have not been characterized. We estimated that the total copy number of the rDNA units in the wheat genome is 11 160, of which 30.5%, 60.9% and 8.6% are located on Nor-B1 (1B), Nor-B2 (6B) and other NORs, respectively. The total length of the NORs is estimated to be 100 Mb, corresponding to approximately 10% of the unassembled portion of the genome not represented in RefSeq v1.0. Four subtypes (S1-S4) of the rDNA units were identified based on differences within the 3' external transcribed spacer regions in Nor-B1 and Nor-B2, and quantitative PCR indicated locus-specific variation in rDNA subtype contents. Expression analyses of rDNA subtypes revealed that S1 was predominantly expressed and S2 weakly expressed, in contrast to the relative abundance of rDNA subtypes in the wheat genome. These results suggest a regulation mechanism of differential rDNA expression based on sequence differences. S3 expression increased in the ditelosomic lines Dt1BL and Dt6BL, suggesting that S3 is subjected to chromosome-mediated silencing. Structural differences were detected in the regions surrounding the NOR among homoeologous chromosomes of groups 1 and 6. The adjacent regions distal to the major NORs were expanded compared with their homoeologous counterparts, and the gene density of these expanded regions was relatively low. We provide evidence that these regions are likely to be important for autoregulation of the associated major NORs as well as silencing of minor NORs.

An efficient approach for the development of genome-specific markers in allohexaploid wheat (Triticum aestivum L.) and its application in the construction of high-density linkage maps of the D genome.

In common wheat, the development of genotyping platforms has been hampered by the large size of the genome, its highly repetitive elements and its allohexaploid nature. However, recent advances in sequencing technology provide opportunities to resolve these difficulties. Using next-generation sequencing and gene-targeting sequence capture, 12,551 nucleotide polymorphisms were detected in the common wheat varieties 'Hatsumochi' and 'Kitahonami' and were assigned to chromosome arms using International Wheat Genome Sequencing Consortium survey sequences. Because the number of markers for D genome chromosomes in commercially available wheat single nucleotide polymorphism arrays is insufficient, we developed markers using a genome-specific amplicon sequencing strategy. Approximately 80% of the designed primers successfully amplified D genome-specific products, suggesting that by concentrating on a specific subgenome, we were able to design successful markers as efficiently as could be done in a diploid species. The newly developed markers were uniformly distributed across the D genome and greatly extended the total coverage. Polymorphisms were surveyed in six varieties, and 31,542 polymorphic sites and 5,986 potential marker sites were detected in the D genome. The marker development and genotyping strategies are cost effective, robust and flexible and may enhance multi-sample studies in the post-genomic era in wheat.

Elucidation of the origin of 'agriocrithon' based on domestication genes questions the hypothesis that Tibet is one of the centers of barley domestication.

Wild barley forms a two-rowed spike with a brittle rachis whereas domesticated barley has two- or six-rowed spikes with a tough rachis. Like domesticated barley, 'agriocrithon' forms a six-rowed spike; however, the spike is brittle as in wild barley, which makes the origin of agriocrithon obscure. Haplotype analysis of the Six-rowed spike 1 (vrs1) and Non-brittle rachis 1 (btr1) and 2 (btr2) genes was conducted to infer the origin of agriocrithon barley. Some agriocrithon barley accessions (eu-agriocrithon) carried Btr1 and Btr2 haplotypes that are not found in any cultivars, implying that they are directly derived from wild barley through a mutation at the vrs1 locus. Other agriocrithon barley accessions (pseudo-agriocrithon) carried Btr1 or Btr2 from cultivated barley, thus implying that they originated from hybridization between six-rowed landraces carrying btr1Btr2 and Btr1btr2 genotypes followed by recombination to produce Btr1Btr2. All materials we collected from Tibet belong to pseudo-agriocrithon and thus do not support the Tibetan Plateau as being a center of barley domestication. Tracing the evolutionary history of these allelic variants revealed that eu-agriocrithon represents six-rowed barley lineages that were selected by early farmers, once in south-eastern Turkmenistan (vrs1.a1) and again in the eastern part of Uzbekistan (vrs1.a4).

Correction to: Generation of expressed sequence tags for discovery of genes responsible for floral traits of Chrysanthemum morifolium by next-generation sequencing technology.

After publication of the original article [1] the authors noted that the following errors had occurred.

Genome-wide analysis of rice cis-natural antisense transcription under cadmium exposure using strand-specific RNA-Seq.

BACKGROUND: The elucidation of novel transcripts and their expression in response to various stress conditions is necessary to understand the transcriptional network of plants as an adaptation to biotic and abiotic stresses. We performed strand-specific RNA-Seq (ssRNA-Seq) on rice exposed to cadmium (Cd) for 24 h and investigated the expression of cis-natural antisense transcripts (cis-NATs), a class of endogenous coding or non-protein-coding RNAs with sequence complementarity to the opposite strands of RAP transcripts. RESULTS: Many RAP transcripts possessed cis-NATs and these cis-NATs were responsive to some extent. Cis-NATs were upregulated from 26, 266 and 409 RAP gene loci, while 2054, 2501 and 2825 RAP transcripts were upregulated from 38,123 RAP loci under high Cd exposure in roots at 1, 12 and 24 h, respectively. In addition, most of the upregulated cis-NATs showed little upregulation under ABA or cold treatment. A number of cis-NATs were upregulated from less than 35 RAP gene loci in different tissue and time-point combinations under low Cd exposure, suggesting that cis-NATs respond to environmental stress. Furthermore, 409 RAP transcripts with upregulated cis-NATs were classified into three groups based on the expression of the RAP transcripts from the opposite DNA strand, including 138 upregulated, 128 invariable, and 143 downregulated transcripts, although the responses of cis-NATs and RAP transcripts were not always correlated. CONCLUSIONS: We have shown that the cis-NATs identified by ssRNA-Seq analysis are novel genes and that some of them are stress-specific and show different responses depending on the degree of stress and tissue. These results improve our understanding of the complete molecular mechanism of plant adaptation to Cd exposure.